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Abstract. In this study we investigate the line-profile variability (Ipv) of Ha for a large sample of 0-type supergiants 
. (15 objects between 04 and 09.7), in an objecti ve, statistical ly rigorous manner. We employed the Temporal 

Variance Spectrum (TVS) analysis, developed by iFuUerton. G ics & Bolton ( 1996) for the case of photospheric 
\ absorption lines and modified by us to take into account the effects of wind emission. By means of a comparative 

■ analysis we were able to put a number of constraints on the properties of this variability - quantified in terms 

^ ' of a mean and a newly defined fractional amplitude of deviations - as a function of stellar and wind parameters. 

^ , The results of our analysis show that all the stars in the sample show evidence of significant Ipv in Ha , mostly 

• dominated by processes in the wind. The variations occur between zero and 0.3 v^o (i.e., below ~ 1.5 i?* ), in good 

I ' agreement with the results from similar studies. 

, A comparison between the observations and corresponding line-profile simulations indicates that for stars with 

1/^ ' intermediate wind densities the properties of the Ha variability can be explained by simple models, consisting 

I of coherent or broken shells (blobs) uniformly distributed over the wind volume, with an intrinsic scatter in the 

. maximum density contrast of about a factor of two. For stars at lower and higher wind densities, on the other 

' hand, we found certain inconsistencies between the observations and our predictions, most importantly concerning 

fH ^ the mean amplitude and the symmetry properties of the TVS. This disagreement might be explained with the 

^ presence of coherent large-scale structures (e.g., CIRs), partly confined in a volume close to the star. 

Interpreted in terms of a variable mass-loss rate, the observed variations of Ha indicate changes of ±4% with 
respect to the mean value of M for stars with stronger winds and of ± 16% for stars with weaker winds. The effect 
of these variations on the corresponding wind momenta is rather insignificant (less than 0.16 dex), increasing only 
^ I the local scatter without affecting the main concept of the Wind Momentum Luminosity Relationship. 

^ ' Key words, stars: early type - stars: mass loss - stars: winds, outflows - stars: activity - methods: data analysis 

. 5t 1- Introduction The most common approach used to study wind vari- 

ability in optical and UV domains is to follow line-profile 
The basic philosophy underlying present day hot star variability (Ipv) of one or several spectral lines, formed 
model atmospheres contains the assumption of a globally in different regions, in order to determine relevant time- 
stationary and spherically symmetric stellar wind with a scales and variability patterns and thus to obtain some 
smooth density stratification. Although these models are insight into the nature and the physical origin of the vari- 
generally quite successful in describing the overall wind ations. This kind of surveys requires long sets of stellar 
properties, there are theoretical considerations, supported spectra with high S/N ratio and high temporal resolu- 
by numerous observational evidences, which indicate that tion, which implies that only few objects have been in- 
hot stars winds are very far from being smooth and sta- vestigated so far. Through such investigations clear evi- 
tionary. dence for the presence of large-scale time-dependent wind 

perturbations (e.g., in the form of Discret e Absorption 

~Z ' ~~ ~~ AT Tvi , Components, DACs) was found in UV fPrinia et al 
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1999; Prima ct al. 2UU2) and optical (I'lilicrton ct ai. 1992; 
Prinia fc F u llcrton 1994: Rauw et al. .2001; Prinia et aL 
200lt lMarkovaH2002r spectra of many O and early B 
stars. Since DACs h ave been observe d in WR stars 
l|Prinia fc SmithI Il992l) and in an LBV llMarkovalllOsl 
as well, they are thought to be a fundamental property of 
radiative driven stellar winds. 

Another source of Ipv in hot stars winds are 
small-scale structures (clumps) which are believed 
to r esult from strong instabi l ities i n the wind it- 
self llOwocki. Ca,stor fc RvbickH EUl iFeldmeieH 
lOwocki fc Puyll999h . While the clumped nature of WR 
winds was unambiguously proven by observations, the 
presence of clumps in 0-star winds has so far re- 
lied on indirect evidence only llCrqwther et al.1 120021: 
iBianchi fc Garciall20n2l:[Markova et al.ll2004^ . 

Wind structures and temporal variability are among 
the most important physical processes that may sig- 
nificantly modify the mass loss rates derived from 
observations. Since accurate mass loss rates are cru- 
cial for evolutionary studies (e.g., iMevnet et alJ Il993) 
and for extra-galactic distance determinations (via the 
Wind Momentum Luminosity Relationship ( WLR) , cf . 
iKudritzki fc Pulj l2000() . it is particularly important to 
know to what extent the outcome of these studies might 
be influenced by uncertainties in Mdu e to the effects o f 
wind structures and variability. Indeed, IKudritzki l|l999|) 
has noted that wind variability is not expected to affect 
the concept of the WLR significantly. However, this sug- 
gestion is based on results obtained via a detailed investi- 
gation of one object alone, while similar data for a large 
number of stars of different spectral types and luminosity 
classes are needed to resolve the problem adequately. 

Following the outlined reasoning, a project to study 
the effects of wind structure and variability in Galactic O- 
type stars has been recentlv started bv our group. While 
in a previous paper ijMarkova et al. ''2004. Paper I) we have 
dealt with problems concerning the WLR and the effects of 
wind clumping, in the present one we address the question 
of wind variability as traced by Ha and the dependence (if 
any) of the properties of this variability on fundamental 
stellar and wind parameters. In particular, in Sect. El we 
describe the observational material and its reduction. In 
Sect. 121 we outline the method used to detect and analyze 
the Ha Ipv. In Sect. ^ Eland the results of our analy- 
sis are presented in detail while in Sect. (HI the outcomes 
of some simple 1- and 2-D simulations are described. In 
Sect. (HI we summarize the major results and give some 
comments and conclusions. 



2. Observations and data reduction 

Our sample consists of 15 Galactic supergiants with spec- 
tral classes from 04 to 09.7, all drawn from the list of 
stars analyzed bv iMarkova et a l. (2004) in terms of their 
mass-loss and wind momentum rates. Table ^ lists the 
objects along with some of their stellar and wind parame- 



ters, as used m the present study. All data are taken from 
Paper I. 

A total of 82 high-quality Ha spectra {R — 
15 000) of the sample stars were collected between 
1997 and 1999. The observations were obtained at 
the Coude focus of the 2m RCC telescope at 
the National Astronomical Observatory (Bulgaria) us- 
ing an ELECTRON CCD (520x580,22x24^) and a 
PHOTOMETRIC CCD (1024xl024,24/x).i For all stars 
but one the S/N ratio, averaged within each spectral 
time series, lies between 150 to 250, while in the case of 
HD 190429 it is -100. 

The temporal sampling of the data for each target is 
not systematic but random, with typical values of the 
minimum and maximum time intervals between succes- 
sive spectra of 1 to 2 and 7 to 8 months, respectively. In 
several cases observations with a time-resolution of 1 to 5 
days are also available, but in none of these cases these ob- 
servations dominate the corresponding time series. Thus, 
we expect the results of our survey to be sensitive to vari- 
ations which occur on a time-scale which is significantly 
larger than the corresponding wind flow time (of the order 
of a couple of hours). 

More information about the ob servational material 
and its reduction can be found in iMarkova fc ValchevI 
(2000) and in Paper I. In particular, to reduce the obser- 
vations we have followed a standard procedure (developed 
in IDL) which includes: bias subtraction, flat-fielding, cos- 
mic ray hits removal, wavelength calibration, correction 
for heliocentric radial velocity, water vapor lines removal 
and re-binning to a step size of 0.2 A per pixel. 

3. Methodology and measurements 

Since we were going to study a large number of objects 
and since in many cases our observations were not sys- 
tematic but with large temporal gaps in between, from 
the onset of this investigation on we recognized that our 
ability to characterize the wind variability of individual 
targets would be restricted, e.g., we would not be able to 
determine time-scales and variability patterns. Moreover, 
to work effectively, we would need to employ some simple 
and fast method both to detect and quantify Ipv and to 
constrain the properties of this variability as a function 
of fundamental stellar and wind parameters of the sample 
sta rs. 

iFullerton. Gies fc BoltonI l)l996l) have developed a very 
sensitive and rigorous method to handle Ipv that is poorly 
sampled in time. This method, called Temporal Variance 
Spectrum (TVS) analysis, was proven to be a powerful 
tool to detect and compare Ipv of pure absorption profiles. 

^ The use of different detectors is not expected to bias the 
homogeneity of our sample because the noise characteristics 
of these two devices are practically the same. The root-mean- 
square (rms) read-out noise of the ELECTRON CCD is 3 elec- 
trons per pixel (i.e 1.5 ADU with 2 electrons per ADU) while 
the rms read-out noise of the PHOTOMETRIC CCD is 3.3 
electrons per pixel (2.7 ADU with 1.21 electrons per ADU). 



lable 1. btcliar and wind parameters oi the sample stars used m the present study. All data are taken Irom 
iMarkova et al. (,2004) . 
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However, its potential with respect to profiles influenced 
by wind emission has not been tested systematically so far. 
Although the application of the TVS technique to detect 
Ipv in emi ssion lines does not seem to pose serious prob- 
lems (e.g- lKaufer et a,].lll"993lKa,ner et a] .111 9971: iMarkoval 
120021) , its implication for the objectives of a comparative 
analysis (e.g. to compare the strength of Ipv in different 
lines or different stars) certainly needs to be carefully in- 
vestigated. 

To study the Ha variability of the stars in our sample 
we modified the main philosophy of the TVS analysis in 
order to take into account the effect of wind emission. 

To compute the TVS of Ha as a function of veloc- 
ity across the line and to determine the velocity width 
over which significant var i ability occurs, AV^ we followed 
iFullerton. Gies fc BoltonI l|l99fit) but assumed that the 
noise is dominated by photon noise. ^ In this case, the TVS 
for the pixels in column j (i.e., at wavelength/velocity j) 
is calculated from 



N 



TVS,=Y. 



w{i){S^j - Sjf 
S,,{N~l) 



(1) 



where Sj is the weighted mean spectrum for the j-th pixel, 
averaged over a time series of N spectra, and given by 



J TV 
with the weighting factors, w^, given by 

= (— ) 



where 



N 



(2) 



(3) 



(4) 



This assumption seems to be justified because we rely on 
Coude spectra of relatively high quality. 



and (Tic is the value of the noise in the i-th spectrum, 
averaged over a certain number of continuum pixels (40 
in our case). 

The RMS deviations {RMS = TVS'^-'^) as a function 
of velocity across Ha for the time series of each target are 
shown in the top panels of Figures ^ to O The level of de- 
viations in the continuum, cto: is represented by a dashed 
line, while the threshold of significant Ipv, fixed at the 
corresponding 99% confidence level of the ctqXn-i distri- 
bution, is marked with a dashed-dotted line. We want to 
stress here that although our implementation is in terms of 
TVS°-^, hereafter we shall continue to refer to the "TVS" 
and the "TVS analysis" , respectively.^ 

To localize the Ha Ipv in velocity space we used the 
"blue" and "red" veloci ty limits of significant var i ability , 
Wb and Vr, introduced bv lPullerton. Gies fc BoltonI l)l99()|) . 
The measurements have been performed interactively to 
fix the positions of the two points where the TVS crosses 
the horizontal line representing the threshold of significant 
Ipv. The accuracy of these measurements depends on the 
quality of the data used and on the strength of Ipv. For 
example, in the limiting case of a strong Ipv (i.e., a TVS 
with large amplitudes and steep spectral gradients) the ac- 
curacy of the individual measurements might be as good 
as ±20 km s~^ . Alternatively, in the case of a weak Ipv 
(e.g., with amplitudes just above the threshold of signif- 
icant variability) the determination of the velocity limits 
might become so uncertain that different positions of al- 
most similar probability may exist for each limit. In these 
latter cases and in order to assess the effects of such un- 
certainties on the outcomes of our analysis, we provide 



Root mean square deviations have been used instead of 
the TVS itself since the former quantity scales linearly with 
the size of the deviations. Thus, it is more appropriate for a 
direct comparison of the strength of Ipv in various stars (see 
also Fullerton et al.). 



two couples 01 estimates tor Vb and v^. itiese two sets oi 
values, expressed in km s^^ , are listed in Column 6 of 
Tablc[51as a first and a second entry. We consider the first 
entry as the more reliable one and will refer to it as the 
"conservative case" . 

As a by-product of the measurement of vi^ and v^, we 
obtain the total velocity width over which significant vari- 
ability in Ha occurs, AV = {vr — Vb). In order to find 
some constraints on the distribution of the Ipv in physi- 
cal space, we furthermore determined the radial distance 
''max, where v{r = Tmax) = Wb, assuming that the wind 
velocity obeys a standard law of the form 



v{r) =Voo{i- , 



1 



Woo 



1//3 



(5) 



(6) 



with P and Vqo from Tableland Umin = 1-0 km s^^ . The 
obtained estimates of r„iax, expressed in units of R^, , are 
given in Column 7 of Table El We are aware of the fact 
that emission variability is difficult to localize and could 
in principle be due to the net effect of fluctuations that 
occur in different locations under different conditions and 
therefore consider these estimates as upper limits only. 

To quantify and compare Ipv, 

iFullerton. Cies fc BoltonI l|l99(t l have introduced two 
parameters, called mean and fractional amplitude of 
deviations, Aip^ and aipv In the following we will refer 
to these quantities as to Ap and op (with "F" referring 
to FuUerton). The first parameter is expressed in units 
of the normalized continuum flux, while the second one 
is a dimensionless quantity. The authors define these 
quantities as follows: 



1 

AV 



TVSf^ dv 



100 /J" {TVSj 



P'' 5,- — 1 dv 



(7) 



(8) 



The above expressions imply the following. First, the 
mean amplitude of deviations is independent of profile 
type and thus can be used for evaluating the statistical 
significance of Ipv both in absorption and in emission pro- 
files. Second, the fractional amplitude depends (via the 
denominator) on the strength of the underlying spectral 
feature but does not make any difference between profiles 
in absorption and in emission. Particularly, it becomes a 
non-monotonic function of wind strength, with a maxi- 
mum in those regions where the wind-emission has (more 
or less) completely filled in the photospheric absorption, 
i.e., where the net equivalent width within [vh, Vr] is close 
to zero. Therefore, this quantity is inappropriate for inves- 
tigating profiles which are influenced by wind emission of 
different extent. Actually, this problem has already been 
outlined by FuUerton et al. 

In order to optimize the fractional amplitude to ac- 
count for the systematic difference in the strength of 



Ha as a function of wind strength, we decided to nor- 
malize the integral over the TVS to a quantity which we 
called "Fractional Emission Equivalent Width" (FEEW). 
With this new definition of the fractional amplitude, now 
denoted by on to distinguish it from FuUerton's parame- 
ter ap, this quantity is a measure of the observed degree of 
variability per unit fractional wind emission. "Fractional" 
refers here to the observed range of significant variability, 
[vii,Vr]. Formally, on is given by'* 



100£{TVS,-aif'dv 



a-N 



phot 



l)dv 



(9) 



The first term in the denominator of EqEI represents the 
fractional equivalent width of the observed profile (posi- 
tive for emission and negative for absorption), while the 
second one gives the fractional equivalent width of the 
photospheric component of Ha (always negative). In total, 
the denominator thus gives the fractional wind emission 
(always positive). A further discussion of an is given in 
Sect. El 

The mean amplitude, as defined by Eq.Q on the other 
hand, does not seem to pose any problem concerning an 
assessement of the statistical significance of variability 
across Ha. In their original study, FuUerton et al. have 
noted that this quantity cannot serve as a comparative 
tool because it does not account for differences in the 
strength of the underlaying absorption feature. In con- 
trast, in the case of Ha from 0-type supergiants the mean 
amplitude might depend on the wind strength^ and might 
therefore become of interest as well, in order to examine 
and to compare the wind variability in stars of various 
spectral types. Motivated by this possibility we re-defined 
the mean amplitude to account (partially) for differences 
in the overall quality (i.e., in S/N) of the time series of the 
sample stars, by subtracting ctq from the TVS, 



An = 



100 
AV 



{TVS, - a'S'' dv. 



(10) 



The photospheric profiles of Ha required to derive the val- 
ues of aN have been selected from a grid of plane-parallel 
models in dependence of the particular stellar parameters 
(Table 1, see also Paper I). Note that the (relative) uncer- 
tainty of the denominator becomes rather large in those 
cases where the wind-emission is only marginal, since in 
this case the errors introduced by uncertainties in the stel- 
lar parameters (affecting the actual choice of the photo- 
spheric profiles) become significant. 

To estimate the uncertai nty in on we followed 
iFullerton. Gies fc BoltonI l)l996|) but used a re-formulation 
(derived by A. FuUerton, priv. com.) of their Eq. 16. 

* In this expression, we have accounted only for uncertainties 
caused by photon noise while the error due to small differences 
in the continuum level of individual spectra in a given time 
series is neglected. 

^ E.g., the numerator in Eq. |71is expected to react on wind 
density (the higher the density the larger the emitting volume) . 



Additionaliy, we assumed that the errors m both ctq and 
FEEW are neghgible and that the accuracy of the devi- 
ations for each pixel j within Ha is identical and equals 
(Tq-^ Under these circumstances, standard error propaga- 
tion gives: 



cr(ajv) 



100 Aw 
FEEW 



1 



2{N ~1) 

0.5 



0.5 



TVSj 



(11) 



where j runs over all the pixels between v\, and fr, while 
-/V and Aw denote the number of spectra in the time series 
and the discretized integration step, respectively. In our 
case Aw ^ 9 km s~^ (not to be confused with the total 
velocity width AV^!). 

The factor of 100 appearing in Eas.lSltollllconverts the 
corresponding quantities to a percentage. The estimates of 
(To, A^r and An ± cr(aN)) for each sample star are listed in 
TableEl Columns 5, 13 and 14, respectively. 

To assess the contribution of changes in Ha line 
strength to the Ipv detected by the TVS analysis, we es- 
timated the mean value of the net wind emission, Wom, 
by subtracting the equivalent width (EW) of the photo- 
spheric profile, tVphot, from the EW of the time-averaged 
observed profiles. 

The observed equivalent widths were measured by in- 
tegrating the line fiux between limits which were set in- 
teractively, judging by eye the extension of the emis- 
sion/absorption wings. These limits did not change for a 
given star but could vary for different stars. The internal 
precision of individ ual EW measurements, est imated in 
the way described in lMarkova & ValchevI l|200flh . is better 
than 10%. The EW of the photospheric component was 
calculated by integrating over the appropriate synthetic 
profile. The Wcm estimates and their (standard) error are 
given m Column 8 of Tab.El together with the EW of the 
photospheric components. 

Since in 0-type supergiants Ha originates from pro- 
cesses taking place in the wind and in the photosphere, 
contributions from absorption Ipv to the observed Ipv 
might be expected (via the photospheric components of 
Ha and Heii A6560 ) . To investigate this possibility we 
consulted the literature (particular references are given 
below) concerning the presence of absorption Ipv in our 
sample. In addition and as a secondary criterion, we used 
the TVS of the Hen A6527 absorption line located at about 
1650 km s~^ blue- wards of Ha. In those cases where the re- 
sults of our TVS analysis of Heii A6527did not agree with 
the results from the literature, the latter were adopted.^ 



The latter assumption is justified since we rely on Coude 
spectra of relatively high signal to noise ratio. 

''' Such inconsistencies may occur because the temporal sam- 
pling of our observations is not well-suited for studying Ipv on 
short time scales (e.g., hours) which seem to be typical for a b- 
sorption Ipv in 0-type stars (|Fullerton. Gies fc Boltonlll99(Jl . 



io obtam constraints on the variability ot mass-ioss 
rates, for each star we determined lower and upper limits 
for M . This has been done by fitting those Ha profiles 
which display the smallest and the largest wind emission 
present in the given time-series, by means of synthetic 
profiles. These have been calculated using stellar and wind 
parameters from Tableland employing the same method 
as used in Paper I. The accuracy of the M determinations 
equals ±20% for stars with H a in emission and ±30 % for 
stars with Ha in absorption l|Markova et al.ll2004l) . The 
estimates of M ^in and M max as well as the amplitude of 
the M variability (given in percent of Mniin) are listed in 
Columns 9 to 11 of Table HJ 

To quantify the wind strength we finally calculated the 
"mean wind density", < p >, using data given in Tabled 
and Table |2|(Column 12), by means of 



<P>-- 



M 



47r(1.4i?0'«o 



(12) 



i.e., we considered the density at a typical location of 
1.4 . 

4. TVS analysis of Ha line-profile variability 

4.1. Stars of early spectral types 

Within our sample, the 04/5 supergiants HD 190429A, 
HD 16 691 and HD 14947 constitute the group with the 
highest effective temperature. Their Ha profiles, displayed 
in FigH appear to be completely in emission and consist 
of a well-developed emission core superimposed onto rel- 
atively strong and extended emission wings. The profiles 
are slightly asymmetric with a red wing being steeper than 
the blue one and a peak emission red-shifted with respect 
to the stellar rest frame. 

For all stars, the mean amplitudes of deviations are 
significant at the 99% confidence level, indicating genuine 
Ipv in Ha. Whereas for HD 190429A and HD 16 691 the 
variations in Ha are distributed preferentially on the blue 
side of the emission peak, for HD 14 947 they extend al- 
most symmetrically with respect to it. For HD 190429A 
two values for Vr are provided, because of the rather small 
amplitudes of deviations at the red edge of the TVS. 

Since the TVS analysis does not show evidence for sig- 
nificant Ipv in Hell A6527(see upper panels of Fig.P), we 
suggest that the variations observed in Ha are mostly (if 
not completely) due to processes in the wind. In the partic- 
ular case of HD 190 429A, this assumption is suppor ted by 
results reported bv lFuUerton. Gies fc Boltm] lll996l) . indi- 
cating that the photospheric ClvA5801 line of this star 
does not show signatures of noticeable Ipv. 

Our measurements indicate that in each of the stars 
considered here the observed Ipv in Ha has been accom- 
panied by real (i.e., extending the measurement errors of 
10%) variations in the equivalent width. In terms of the 
adopted model these variations in the Ha line strength 
can be reproduced by variations in M ranging from 8% 
(HD 190429A) to 10% (HD 14947). 
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Fig. 1. Stars of early spectral type (04/5). Upper part of each panel: Mean Ha profiles (thick line) and RM S deviations 
as a function of velocity across the line. Middle part of each panel: Time-series of observed Ha profiles. Lower part of 
each panel: Ha profiles with maximum and minimum wind emission in the time series. Velocity scale centered at the 
corresponding systemic velocity. 



4.2. Stars of intermediate spectral type 

The subset of our supergiant sample with intermediate 
effective temperatures (32 000 K < Tcff < 37000 K) in- 
cludes six objects with spectral types ranging from 06 to 
08.5. Apart from HD 24912, whose Ha is in absorption 
and which will be considered separately at the end of this 
subsection, all the other objects exhibit Ha in emission 
(see Fig. El). 

The observations indicate the presence of noticeable 
Ha Ipv in any star of this sub-group (see the middle parts 
of the panels in Fig. EJ. The emission component varies 
from sharp, single-peaked emission to a double-peaked and 
somewhat broader structure, while the absorption compo- 
nent changes from a well-defined trough to a weak, hardly 
noticeable feature. At the same time, the morphology of 
the profiles does not change drastically and always consists 
of a P Cygni-like core on top of extended emission wings. 
The emission peak of the mean Ha profiles is red-shifted 
with respect to Ugys , whereas the position of the absorp- 
tion dip is different for various stars and ranges from 260 
to 360 km s~^ relative to the emission peak. 

For the majority of the stars our TVS analysis re- 
veals the presence of significant Ipv concentrated at the 
P Cygni-like core. The variations are distributed almost 
symmetrically with respect to the rest wavelength. Only 
in HD 225 160 and in HD 338 926 they are concentrated 
towards the blue and the red side of the line, respectively. 
The velocity limits of significant variability are always well 
defined except for HD 210 839, where these limits seemed 
to be more uncertain. 

To clarify the situation, we investigated the individual 
Ha profiles of the HD 210839 time-series and found that 
part of the Ipv detected between -900 and -500 km s~^ and 



between 500 and 700 km s^^ might be due to imperfect 
telluric line correction. This possibility was taken into ac- 
count by providing two entries for Wb and Vi- of HD 210 839 
in Table 13 

Interestingly, few years ago iKaoer et alJ l|l999l) re- 
ported evidence for significant variability in Ha of 
HD 210839 concentrated towards the blue side of the line 
(between and -400 km s^^ ). The authors attributed 
this variability to the cyclic (P = 1.4 days) appearance 
of DACs in UV resonance lines. This finding is somewhat 
different from the one derived by us, where significant vari- 
ability in Ha of HD 210 839 extends towards the emission 
peak of the line as well. This "inconsistency" might be ex- 
plained by suggesting that in addition to the short-term 
variability caused by DACs there is another variability 
component that operates on a time scale longer than cov- 
ered by the observations used by Kaper et al. 

Interestingly as well, the morphology of Ha of 
HD 192 639 does not indicate any evidence for deviations 
from spherical symmetry and hom ogeneity, in contrast to 
the findings bv lRauw et alJ l)200l|) . If not due to observa- 
tional selection, this r esult might s ugges t that the wind 
structure observed by iR.aiiw et al.l l|200lh was not a per- 
manent but a transient feature of the wind, similar to the 
one observed in a Cam (Markov a 2002). 

Except for HD 17603 and HD 338 926 no indica- 
tion of significant absorption Ipv was found. Regarding 
HD 210 839 and HD 19 2 639 this finding does no t agree 
with re sults reported bv iFiillerton. Gies fc BoltonI ill 996h 
and by Ide Jong et al.l l)l999() .^ Thus, we conclude that at 
least for four out of the five stars in this group contribu- 
tions from photospheric Ipv can be expected. 

* Ide Jong et al.l 1^9^ have recognized HD 210 839 as a non- 
radial pulsator. 
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Fig. 2. As Fig.n but for stars of intermediate spectral type (06 to 08.5). 



For all stars a genuine variability in the equivalent 
width of Ha has been found. Neglecting the possible con- 
tribution from absorption Ipv (if present) we found that 
the (rather extreme) variations in the Ha net wind emis- 
sion can be accounted for by a 12% (HD 225 160) to 44% 
(HD 210839) variation in M . 

HD 24 912 Per) was observed once in 1997, three times 
in 1998 and once in 1999. The Ha profiles (lower-right 
panel of Fig. [3) appear completely in absorption, in con- 
trast to the profiles of the other sample stars of similar 
spectral type. The red wing is steeper than the blue one 
and the absorption dip is blue-shifted with respect to the 
stellar rest frame. Profiles with similar signatures are typi- 
cal for stars with weaker winds where the Ha photospheric 
absorption is partly filled in by wind emission. 

The TVS of Ha consists of a blue-shifted, single- 
peaked component plus a double-peaked structure with 
maximum amplitudes concentrated at the absorption core 



and at the red extension of the profile. In an absorp- 
tion profile, a double-peaked TVS might indicate the 
presen ce of radial velocity variabilit y caused by pulsa- 
tions l|Fullerton. Oies fc Boltonlll996l) . Indeed, our mea- 
surements show that the velocity of the absorption core 
of Ha varies around a mean value of -43±13km s~^^. 
Since the position of the absorption dip does not seem 
to depend on the line strength (i.e., on the strength of 
the wind emission) we suggest that the observed radial 
velocity variability is mostly (if not completely) due to 
changes in th e stellar phot osphere and probably caused 
by pulsations l|de Jong et al... 1999. .2001) . 

In addition to this, another variability component 
seems to be present in HD 24912, as indicated by the 
blue-shifted, single-peaked feature in the Ha TVS, located 
between -200 and -450 km s^^ . This velocity interval is 
similar to the interval of the 2-d period variation estab- 

® The uncertainty in individual radial velocity measurements 
equals to ±4.5 km s~^ , i.e., half the bin step of 0.2 A. 
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Fig. 3. As Fig.m21 but for stars of later spectral type (09 to 09.7). Note the single-peaked, blue-shifted absorption 
pr ofiles observed in HP 188209 in June, 1997 and February, 1998. Profiles with similar shape have not been observed 
bv llsraelian et all l|200(/) throughout their long-term monitoring campaign. Note also the P Cygni-like profile observed 
in HD 18 409 in June, 1998. This profile is completely different from the rest of the time series and indicates a strong 
increase in density in the innermost wind region. 



fished bv lKaper etlll l)l997|) and bv lde Jong et all l)200l|) . 

and both findings might have the same origin. 

The EW of Ha varies between 1.3 an d 1.8 A, in 
good agreement with the limits derived by iKaper et alJ 
l|1997(l . Interpreted as due to variations in the mean wind 
density, these limits comply with a ±16% variation in 
Af , which is smalle r than the error of the individual 
M estimates (±25%. 'Markov a et al.ll2004 and thus in- 
significant. Again, this find ing is consi stent with corre- 
sponding results from lkaper et alJ (|l993). Based on simul- 
taneous UV and optical observations of ^ Per, the latter 
authors suggested that the EW variability of Ha is caused 
by the presence of large-scale, time-dependent structures 
in the wind. 



In summary, we conclude that the variability we have 
observed in ^ Per is a mixture of variations originating 
both from the photosphere (caused by non-radial pulsa- 
tions) and from the wind (presumably connected to the 
appearance of DACs). 



4.3. Stars of late spectral type 

The subset of sample stars of late spectral type (09 to 
09.7, Tcff < 32 000 K) includes 6 objects. Five of them 
show Ha profiles with similar morphology whereas an- 
other one - HD 210809 - exhibits a multitude of differently 
shaped Ha profiles and will be considered separately at the 
end of this subsection. 



rtie tia prohics ot HU 188 209, BD+56 739, 
HD 209 975, HD 218 915 and HD 18 409 appear in ab- 
sorption, partly filled in by wind emission (see Fig. 
The shape of the profiles can vary, from star to star and 
for a given star as a function of time, from double-peaked 
absorption with a central reversal peaked at the stellar 
rest frame to an asymmetric, blue-shifted absorption fea- 
ture, with the red wing being steeper than the blue one. 
Ha profiles with similar signatures have been fou nd in 
early B-type supergiants as well, e.g.. lEbbet3 

No indication of extended emission wings has been 
found in this subgroup, except for BD +56 739. This ob- 
ject shows weak, but clearly visible emission wings extend- 
ing to about ±1200 km s~^ , suggesting a relatively strong 
wind."'^^ 

All stars show evidence of real Ipv in Ha. For part 
of the sample stars, the deviations are distributed almost 
symmetrically with respect to the stellar rest frame, with 
maximum amplitudes concentrated almost at the central 
reversal, whereas for others (e.g. HD 209 975), the varia- 
tions are stronger blueward from the rest wavelength. In 
the particular case of HD 209 975 this blue-to-red asym- 
metry of the TVS might be expl ained by bluewards mi- 
grating DACs llKaper et al.lll997l) . Two entries are given 
for HD 18409, since the blue velocity limit for significant 
Ha variability is somewhat uncertain. 

For none of the stars in this group we found evidence 
of significant variability in the Hell A6527 absorption 
line . This finding is consisten t with the results reported 
by iF^iUerton. Gies fc BoltonI l|l99fil) for the two stars 
in common, H D 188 209 and HD 209 975. In contrast, 
llsraelian et al.l ifioOO ^ have reported evidence of quasi- 
periodic absorption Ipv (with P = 6.4'^) for HP 188 2 09. 
Since the data-set used by llsraelian et al.l l)2000|) is 
much more extended th a n the one used by us and by 
iFullerton. Gies fc BoltonI l|l99(j) . we consider their result 
as more reliable. Thus we assume that in five out of the 
six stars the observed variability in Ha is dominated by 
changes in the wind and that only in HD 188 209 a con- 
tribution from photospheric Ipv is to be expected. 

Our measurements show that the main source of Ipv 
in Ha are changes in equivalent width. For all stars the 
M variations needed to account for the extreme changes 
detected in Wem exceed the error of individual determina- 
tions and are therefore considered as real. As an example, 
for HD 188 209 we derived upper and lower limits of 1.5 
and 1.75 xIO^^A/q per year, res pectively, in full agree ment 
with the estimates reported bv llsraelian et all l|200nh . 

HD 210809. The individual spectra, shown in the mid- 
dle part of the corresponding plot in Fig. (31 indicate the 
presence of dramatic Ipv in Ha. The profiles change from 
a relatively weak absorption trough with a flat core via 
an ordinary /reverse P Cygni-like feature to a triple emis- 
sion structure. In addition, strong emission wings extend- 



Weak emission wings (±900 km s ^ 
in HD 209 975. 



might be present also 



mg to ^1500 km s^ are clearly visible. Line profile vari- 
ability with similar signatures has been discussed as an 
indication of long-lived, large-scale wind density pertur- 
bation(s), which co-rotate with the star, giving ri se to ad- 
dition al line emission at various frequencies (;I!.auw et alJ 
I2OOII and references therein). 

As might be expected from the observations, the dis- 
tribution of Ipv in Ha is asymmetric with respect to the 
rest wavelength, with maximum amplitudes concentrated 
at the P Cygni-like core. The smaller amplitude devia- 
tions located between 600 and 900 km s^^ are caused by 
the appearance of a bump on the red emission wing in the 
June 1998 line profile. 

From the TVS of Hen A6527 , we found no indication of 
significant photospheric Ipv, suggesting that the observed 
variations in Ha are caused by changes in the wind. 

Although the observations give clear evidence for de- 
viations from spherical symmetry, we applied our line- 
synthesis code (based on a spherical model) to find con- 
straints on the mass-loss rate variability of the star. 
Fitting particularly the first and the last profile of the 
time series by model calculations, we found that the ex- 
treme variations in Wom can be reproduced by variations 
of ±20% in M . 

5. Ha line-profile variability as a function of 
stellar and wind parameters 

In order to obtain further clues concerning the origin of 
wind variability (as traced by Ha) in O supergiants, we 
examined various correlations between line profile param- 
eters and parameters of the TVS of Ha, on the one hand, 
and fundamental stellar and wind parameters of the sam- 
ple stars, on the other. To search for such correlations, we 
used the Spear man rank-order corr elation test, described, 
for example, bv lPress et al The main advantage 

of this test is that in addition to the correlation coeffi- 
cient (more precisely, the linear correlation coefficient of 
ranks) it also calculates the statistical significance of this 
correlation (expressed as the two-sided significance of its 
deviation from zero), without any assumption concerning 
the distribution of uncertainties in the individual quanti- 
ties. 



5.1. Ha profile shape as a function of spectral type 

An inspection of the mean Ha profiles of the sample stars 
(all supergiants!), displayed in Figures ^^00 shows that 
these profiles evolve as a function of spectral type from 
a slightly asymmetric emission with a peak value red- 
shifted with respect to Usys , via an emission feature with 
a P Cygni-like core, to a feature in absorption (with or 
without central emission reversal). In stars of early and 
intermediate spectral type extended emission wings can 
be seen, while in stars of late spectral type the presence 
of such wings is rare. 

There are two stars that deviate from this behaviour, 
HD 24 912 and HD 210809. The former one exhibits a 



pure absorption proiiie instead oi a Cygni-like prohle 
(see Fig. Consequently, its Ha line resembles much 
more those profiles from luminosity class III than fr'om 
luminosity class I objects of the same spectral type. As 
we have already pointed out in Paper I, the parameters 
of HD 24 912 are somewhat insecure due to the uncertain 
distance - H D 24912 is not a member of PerOB2 but a 
runaway star l|Gieslll98^ . Thus it is rather likely that the 
discrepancy in profile shape is due to an erroneous (re- 
) assignment in luminosity class l|Herrero et al ]Ll992) and 
sugges t that the original value as assigned by .Walborni 
lll973l). luminosity cl ass III, is more appropriate (see also 
iRepolust et"al]l2004j) . 

The second outlier, HD 210 809, shows a P Cygni- 
like profile instead of an absorption profile partly filled 
in by wind emission. This is the only star in the sam- 
ple for which our observations suggest a strong devia- 
tion from spherical symmetry. From the similarity to the 
Ha and Heii 46 86 time-serie s of H D 192 639 observed 
and discussed by iRauw et alJ l)200l|) (single and double 
peak structure in emission), we speculate that also here 
a "confined co-rotating wind" is present. This interpreta- 
tion, if correct, would explain the "peculiar" shape of the 
mean Ha profile derived by us. Hereafter, we will refer to 
HD 24912 and HD 210 809 as to "peculiar" stars. 

The observed evolution of Ha in 0-supergiants with 
spectral type (actually with Toff ) is in fair agreement 
with results from theoretical line-profile computations per- 
formed in terms of NLTE, spherically symmetric, smooth 
stellar wind models, although the strength of the observed 
P Cygni-like core can not be reproduced in most cases (e.g., 
iRepolust et"al]l2004() ." The main drivers of this evolu- 
tion are: decreasing line emission caused by decreasing 
wind density (since M decrease s with decreasing Tcs and 
logL/L0 , see lVink et aLlliooilll and decreasing contribu- 
tion of the Hell A6560 blend. (In contrast, purely pho- 
tospheric Ha profiles of a given luminosity class do not 
change significantly as a function of Teff , because in this 
temperature regime the photospheric ionization fraction 
of neutral hydrogen remains fairly constant.) Outliers can 
occur either as a result of strong deviations from spheri- 
cal symmetry and homogeneity in the wind (due to, e.g., 
fast rotation, CIRs, clumps) or as a result of an erro- 
neous spectral type/luminosity class classification or un- 
certain/wrong parameters. 

5.2. Red-shifted emission-peaks 

Our observations suggest that the position of the emission 
peak of the Ha profiles of O supergiants depends on the 
strength of the wind: for stars with weaker winds (Ha in 
absorption with/ without central reversal) this peak is cen- 

This problem has not been solved satisfactorily so far. In 
particular, it is still unclear if the blue-shifted absorption core 
is solely related to the Hell blend or other effects (clumps, 
deviations from 1-D geometry) have to be accounted for addi- 
tionally. 



tered almost at the rest wavelength, whereas for stars with 
stronger winds (Ha in emission) it is red-shifted instead. 
This observation is supported by results of the correlation 
analysis, which shows that the velovity of the emission 
peak, Wo correlates significantly with < p >, (.79/0.0007) 
and in addition with Toff (. 82/0.0003). 

Hereafter, numbers in brackets denote the Spearman 
rank correlation coefficient and the two-sided significance 
of its deviations from zero. Since the latter quantity mea- 
sures the probability to derive a given correlation coeffi- 
cient from uncorrelated data, smaller values mean higher 
significance of the correlation. 

Red-shifted emission peaks have been observed in UV 
resonance lines of 0-type stars, where this finding (in par- 
allel with the presence of an extended absorption trough) 
can been explained in terms of "micro- turbulence" effe cts, 
with Wmicro of the o rder of 0.1 v^o (e.g..lHamannlll980l and 
references therein, iGroenewegen fc Lamer j Il989|) . On a 
first glance, the phenomenon seen in Ha seems to be some- 
what similar. In contrast to the situation for UV resonance 
lines, however, our Ha profile simulations meet no prob- 
lem in reproducing the red-shifted peak, even if the shift 
is large, without any inclusion of micro-turbulence. This 
can be seen clearly by comparing t heoret i cal profiles with 
observ ations, e.g., iMarkova et alJ l|2004|) : IRepolust et alJ 
l)2004l) . A closer inspection of the profile formation pro- 
cess reveals that the apparent shift of the emission peak 
results (at least in our simulations) from the interaction 
between the red-side of the Stark-broadened photospheric 
profile and the wind emission. Let us note that we do not 
exclude the presence of micro-turbulence but that we sim- 
ply do not need it to reproduce the observed amount of 

Wo. 

5.3. Properties of the TVS as a function of stellar and 
wind parameters 

Before investigating the properties of the Ha TVS for our 
sample stars, let us point out that all results outlined be- 
low refer to the "conservative case" (see Sect. 13). Although 
the "non-conservative" data have not been analyzed in de- 
tail, they are included in the corresponding plots and we 
will comment on their infiuence on the final outcome. 

5.3.1. Distribution of Ha line-profile variability in 
velocity space 

In Figure^ we show the blue and red velocity limits (left 
panel, dashed) and the velocity width AV^ (absolute value, 
right panel) of Ha Ipv as a function of the mean wind 
density. Thick vertical lines correspond to the projected 
rotational speed, itusini. Asterisks refer to the conser- 
vative estimates for AV^, while diamonds mark the non- 



In this particular case and because of the reasons outlined 
above, the "pecuhar" stars HD 210 809 and HD 24 912 have 
been discarded from the correlation analysis. 
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Fig. 4. Blue and red velocity limits (left panel) and velocity width (right panel) of significant Ipv in Ha, as a function 
of mean wind density of the sample stars. Thick vertical lines denote the corresponding projected rotational velocity, 
zbwsini. Asterisks refer to the conservative estimates for AV^, while diamonds mark the non-conservative ones. 

Table 2. Ha line profile and variability parameters. N denotes the number of available spectra. Vc is the velocity 
of the emission peak while Vh, Vr are the "blue" and "red" velocity limits of significant variability. All velocity data 
are measured with respect to the stellar rest frame and given in km . ctq is the standardized dispersion of the 
corresponding time-series, rmax (expressed in R^, ) denotes the upper limit in physical space where significant variations 
in Ha are present. Wem is the mean equivalent width of net wind emission and its standard deviation, both given 
in A. Wphot is the equivalent width of the photospheric component. Af and Mmax (in 10~^ Mq /yr) denote the 
corresponding limits if the observed variability is attributed to variations in M alone, while AM is the amplitude of 
this variability expressed in percents of M min- < p > is the mean wind density fEa. ll2|l . and A-^ and cn are the mean 
and the fractional amplitudes (Eas. [TUlandl^. respectively. 
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HD 190 429 
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HD 16 691 
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[ -448, 266] 


1.22 


10.87±0.81/3.23 


12.0 


13.0 


8% 


13.13 
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7.17± 0.02 


HD 14 947 
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4 


no 


82 


0.50 


[ -704, 747] 


1.43 


9.35±0.90/3.24 


14.5 


16.0 


10% 


13.27 


1.90 


7.84± 0.02 


HD 210 839 


4 


11 


yes 


220 


0.67 


[ -633, 478] 


1.40 


4.74±0.55/3.00 


6.8 


9.8 


44% 


13.42 


1.92 


11.63± 0.03 
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1.73 
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HD 192 639 


5 


7 


no 


140 


0.67 


[ -393, 460] 


1.27 


6.15±0.52/3.02 


4.7 


5.4 


16% 


13.38 


3.36 


13.77± 0.03 


HD 17 603 


6 


7 


yes 


210 


0.60 


[ -373, 619] 


1.27 


4.04±0.49/2.84 


5.5 


7.2 


31% 


13.56 


1.87 


12.49± 0.05 


HD 24 912 


7 


5 


yes 


-50 


0.40 


[ -419, 309] 


1.12 


1.36±0.14/2.85 


4.5 


5.2 


16% 


13.77 


1.09 


15.32± 0.14 


HD 225 160 


8 


4 


no 


140 


0.46 


[ -491, 384] 


1.33 


4.51±1.04/2.62 


5.1 


5.7 


12% 


13.45 


3.21 


15.61± 0.02 


HD 338 926 


9 


3 


yes 


140 


0.44 


[ -352, 362] 


1.21 


4.54±0.38/2.61 


4.5 


5.4 


20% 


13.59 


2.48 


11.78± 0.03 


HD 210 809 


10 


5 


no 


10 


0.63 


[ -336, 909] 


1.15 


3.63±0.63/2.59 


3.2 


4.5 


41% 


13.60 


2.77 


24.11± 0.05 


HD 188 209 


11 


6 


yes 


-6 


0.72 


[ -176, 144] 


1.09 


1.56±0.19/2.18 


1.5 
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17% 


13.87 


2.67 
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BD+56 739 


12 


3 
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1.12 
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2.5 


19% 


13.80 
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13 
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no 
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0.48 
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1.09 


1.36±0.29/2.18 


1.5 


1.9 


27% 


13.92 


1.90 


20.85± 0.10 


HD 218 915 


14 


6 


no 


-8 


0.70 


[ -149, 206] 


1.07 


1.81±0.24/2.18 


1.6 


2.0 


25% 


13.91 


2.51 


17.14± 0.10 


HD 18 409 


15 


5 


no 


3 


0.70 


[ -274, 268] 


1.08 


1.40±0.47/2.17 


1.5 


2.2 


47% 


13.64 


3.39 


47.18± 0.11 



conservative cases. In combination with the results de- 
scribed in Sect.0I these figures indicate that: 

i) for all stars the Ha Ipv extends beyond the limits de- 
termined by stellar rotation and thus must be linked 
to the wind; 

ii) in most of our sample stars the variations occur either 
symmetrically (within the error) with respect to the 
rest wavelength or with a weak blue-to-red asymmetry. 
For two objects, HD 17 603 and HD 210809, the TVS 
has a noticeable red-to-blue asymmetry, while in other 



two stars, HD 190429A and HD 16 691, the variations 
are stronger and more extended bluewards of the rest 
wavelength. 

iii) the velocity width for significant variability in Ha is 
larger in stronger winds than in weaker ones. There 
are two stars that deviate from this rule: HD 190429A 
and HD 16 691 which exhibit variations over a veloc- 
ity interval that is considerably smaller than expected 
from the strength of their winds. 
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Fig. 5. Examples of significant correlations between the fractional amplitude of deviations, cn, as defined in the present 
study, and stellar and wind parameters of the sample stars. Distributions of the mean amplitude An are shown for 
comparison. Asterisks refer to the conservative estimates, while diamonds mark the non-conservative ones. Positions 
of the two "peculiar" stars are denoted by 'P'. 



Further analysis of the velocity data listed in Table |21 
shows that in all sample stars significant Ipv in Ha occurs 
below 0.3 Voo ■ Converted to physical space, this yields an 
upper limit of roughly 1.5 i?^for the observed variability. 
Additionally, we found a significant correlation between 
Tmax and log < p >, (0.87/0.00001). This result as well 
as the possible dependence between AV and log < p > 
(0.62±0.01) are readily understood in terms of an increas- 
ing wind volume which contributes to the Ha emission, as 
a function of wind density. 

We are aware of the fact that due to observational 
selection effects and other uncertainties affecting the de- 
termination of the velocity limits, the results described 
above might be questioned. However, note that: (i) the 
probability to obtain a symmetric TVS for a star with a 
strongly asymmetric wind, using snapshot observations, 
is very low; (ii) the established correlations between AV 
and log < p > (right panel of Figure 0)) and between 
'"max and log < p > are both physically reasonable, which 
in turn supports the reliability of the limits determined 
by us; (iii) although the limits and hence AV changes 
when non-conservative instead of conservative estimates 



are considered, the final outcome does not change (see 
Figure Hand Table EJ. 

Thus, we assume that the velocity limits determined 
here do not seem to be strongly biased by either observa- 
tional selection or uncertainties in the measured quanti- 
ties. (But see also the next sub-section.) 

If a good temporal resolution of the variability time- 
scale is provided, the velocity distribution of Ipv in Ha will 
allow to obtain signifi cant information about the wind ge- 
ometry l|Harriesl20nd) . On the other hand, and as we shall 
show later on (see Sect , even in the case of snapshot 
observations some hints about wind structures can be de- 
rived. 

5.3.2. Mean and fractional amplitudes of deviations as 
a function of stellar and wind parameters 

Our TVS analysis shows that the mean amplitude of devi- 
ations always exceeds the corresponding threshold for sig- 
nificant variability, indicating genuine variability in Ha. 
The actual values of An range between 1 and 4 percent of 
the continuum flux, without any clear evidence for depen- 
dence on stellar and wind parameters of the sample stars. 



ihis result is iliustratcd m the upper lett panel ol rig-loj 
where the esthxiates for An are shown vs. log < p >. 

This independence of An on log < p > can have a 
twofold interpretation: First, if the mean amplitude is a 
reasonable measure for wind variability, then the wind 
variability is actually more or less independent on wind- 
strength. Second, the mean amplitude is not the appro- 
priate tool to compare the strengths of Ha Ipv. 

Leaving aside these two possibilities, let us first con- 
sider the following problem. If we assume that the sources 
of observable variability are distributed over a certain vol- 
ume which increases as a function of mean wind den- 
sity (as it is suggested from the increase of rmax with 
log < p >), then one should expect that also An should 
increase with mean density, since the numerator of this 
quantity (the integral over {TVS — cr^)"'^) increases as a 
function of the emitting volume, whereas the denominator 
corresponds to an (increasing) 1-D quantity only. The fact 
that the observed mean amplitude is actually independent 
on log < p > shows that such a simple model is not suf- 
ficient to explain the observations. We will come back to 
this point again in Sect.lHI 

In contrast to the established independence of the 
mean amplitude on wind density, our analysis shows 
the presence of a negative correlation between the frac- 
tional amplitude of deviations, cn, and a number of stel- 
lar/wind parameters. Scatter plots for the strongest cor- 
relations, with Teff (0.92/0.000001) and with log < p > 
(0.80/0.0003), are illustrated in Fig. El In particular, the 
decrease of on with increasing log < p > suggests that 
the observed variability per unit fractional net emission is 
smaller in denser winds than in thinner ones. 

The reliability of the derived values of the mean and 
fractional amplitudes has been checked in two ways: first, 
we examined the stability of the results against effects 
caused by observational selection and other uncertainties 
in the measured quantities; second, we checked the va- 
lidity of the assumptions underlying our definitions of An 
and An. 

In particular, to clarify to what extent these quanti- 
ties (and again AV) might be influenced by observational 
selection effects, we proceeded as follows: 

i) For the two stars with longest time series we reduced 
the number of spectra by about 30% while keeping the 
two most "extreme" profiles. The effect of this data 
manipulation on the TVS parameters was found to be 
insignificant. 

ii) For the stars with longer time series {N > 6) we re- 
moved: (a) the spectrum with minimum Ha emission; 
(b) the spectrum with maximum Ha emission and (c) 
the two spectra with minimum and maximum wind 
emission and re-calculated the TVS. In all the three 
cases the new estimates of ctq and AV turned out to 
be quite similar to the original ones: the established 
differences were less than a few percent. At the same 
time the TVS amplitudes changed by less than 10 per- 
cent for stars with Ha in emission and by about 15 



to zO percent for stars with Ha m (partly re-hiled) 
absorption. The effect of these changes was again in- 
significant concerning the final results for AV, An and 
On- 

Let us finally point out that the estimates of the TVS 
parameters are expected to depend strongly on the qual- 
ity of the data used (i.e., on (To)- To simulate higher 
S/N (about two times higher than the original values) we 
smoothed the spectra in each time series, using a boxcar 
average with a width of 4 pixels, and analyzed them in 
the same way as the original spectral series. Interestingly, 
while the effect of this data manipulation on the estimates 
of AV^ was surprisingly small (less than ±12 percent of 
the original values), the reaction of the TVS amplitudes 
turned out to be quite strong: both An and on, averaged 
over the whole sample, decreased by 37 and 33 percent, 
respectively, compared to their original values. -'^^ Most im- 
portantly, however, the new estimates of AV, An and an 
were found to obey similar dependences on log < p > as 
implied by the original data set.-'^^ 

Summarizing we conclude that the derived indepen- 
dence of An on stellar and wind parameters as well as 
the negative correlations between these parameters and 
On cannot be explained in terms of either observational 
selection or uncertainties in the measured quantities. 

6. Simulations of Ipv in Ha 

In order to account for the systematic difference in the 
strength of Ha as a function of spectral type/mean wind 
density, in Sect.|21we optimized the fractional amplitude 
of deviations by normalizing the integral over the TVS to 
a quantity which we called fractional emission equivalent 
width, FEEW (see Eq.EJ. 

The parameter an defined in this way is thus a mea- 
sure for the observed variability (represented by the cor- 
responding TVS) per unit fractional wind emission and 
has been used in the previous section to investigate the 
dependence of the observed variability on wind density. 
The results obtained might be interpreted as an indica- 
tion that denser winds are less active than thinner winds, 
a finding which would give firm constraints on present hy- 
drodynamical simulations. 

Note, however, that (i) our definition of on implicitly 
assumes that the TVS amplitude is proportional to the 
corresponding amount of wind-emission and that (ii) this 
assumption has not been checked so far. In particular, if 
this assumption was justified, the derived values of cn 
would provide a robust measure for the "observed" degree 
of wind-variability, as it is true for the photospheric Ipv's 
described in terms of ap. 

Part of this decrease might be due to the fact that the 
contribution of higher frequency variability (if any) has been 
reduced by smoothing. 

This result might no longer be true if large differences in 
the quality, i.e., in the standardized dispersion of the spectral 
time-series, ctq, were present (the requirement of homogeneity). 



lable 6. bummary or simple 1-U simulations. Models de- 
noted by"SS" refer to spherical shells, models denoted by 
"BS" to broken shells, respectively. 



Series 


properties 




max((5p//9o) 


SSI 


5w = 0.5wth(H) 




±0.7 


SS2 


5i; = 1.0t;th(H) 




±0.7 


BSll 


Sm = const, Ap(core) 


= 0.1 i?* 


±0.35 


BS12 


Sm = const, Ap(core) 


= 0.1 


±0.7 


BS21 


5m = const, Ap(core) 


= R* 


±0.35 


BS22 


5m = const, Ap(core) 


= R* 


±0.7 


SS3 


5m = const 




±0.35 



6.1. 1-D model simulations 

Thus, a test of our hypothesis is urgently required. Ideally, 
such a test would make use of at least 2-D models of insta- 
ble winds, since the assumption of 1-D shells most proba- 
bly overestimates the actual degree of variability. 

Since 2-D simulations involving a consistent physical 
descr iption are just at their beginning l|Dessart fc Owock3 
l200,Sl) and since, to our knowledge, even for somewhat 
simpler multi-dimensional models no investigation con- 
cerning the dependence on wind parameters is available 
(Sect. 1^?^ . we have proceeded in the following way. 

We have constructed a large number of very simple 
wind-models with variable Ha wind emission, in complete 
analogy to our stationary description (Paper I). The re- 
sulting profiles (10 per model) have been analyzed in the 
same way as the observed ones, i.e., by means of the TVS- 
analysis as described in Sect. 01 To allow for a direct 
comparison with results from our observations, we have 
added artificial Gaussian noise to the synthesized profiles 
(S/N — 200, which is a typical value), and have re-sampled 
the synthetic output onto constant wavelength bins corre- 
sponding to an average resolution of 15 000. 

In order to account for the effects of wind disturbances 
of different size and density contrast in different geome- 
tries, we calculated various series of models: three consist- 
ing of spherical shells (series SS) and four consisting of 
broken shells (i.e., clumps, series BSl and BS2). A sum- 
mary of the various models and their designation is given 
in Table El 

All simulations are based on our (quiet) model for 
HD 188 209 withM = l.e-lQ-SMg/yr (cf. Paper I). In or- 
der to investigate the reaction of the TVS of the synthetic 
profiles as a function of wind-strength we calculated, for 
each model series, 9 different models, with AlogM « 
0.1, particularly at 0.8,1.25,1.6,2.0,2.5,3.2,4.0,5.0 and 
10.0 • lO~^M0/yr. In this way, profile shapes varying from 
pure absorption via P Cygni type to pure emission have 
been obtained, covering all "observed" mean wind densi- 
ties. 

In all these models, only the density was allowed to 
be variable, whereas the velocity field and thus the NLTE 
departure coefficients (as a function of velocity) have been 
kept at their original value. 



In a hrst step (series bbl/bb2), we used the most sim- 
plistic approach of disturbing the density, namely we var- 
ied this quantity in a random way as a function of ra- 
dius only, i.e., we assumed spherical shells. This step, al- 
though rather unrealistic, has been performed particularly 
to check the stability of results from somewhat more "so- 
phisticated" models (series BS1/BS2), which are described 
below. 

The variations are defined in such a way as to allow 
for both positive and negative disturbances around the 
"quiet" model, in order to preserve the mean profile. We 
divided the wind into shells of equidistant velocity range, 
Svsheii, where 

SVsheU=C-Vth(il), c= 0.5,1 (13) 

with Wth(II) the thermal velocity of hydrogen. The differ- 
ent multipliers c define two different series, SSI and SS2, 
respectively. Inside each of the shells, the density has been 
perturbed by a maximum amplitude of ±70%, 

p = Po{l + 6p/po), 5p/po = -0.7+1.4-RAN, (14) 

with pq the stationary density and ran a random number 
uniformly drawn from the interval [0,1]. The specific max- 
imum amplitudes (for series SS, but also for series BS, see 
below) have been chosen in such a way that the resulting 
TVS-integrals and FEEWs are (roughly) consistent with 
the observed values. Lower maximum amplitudes would 
result in a too low degree of variability, and higher ones 
in too large values. 

Note that the density contrast has been assumed to 
be constant within each of the shells, i.e., the number of 
drawn variables is given by the total number of shells, 
which is of the order of 80 for Voo — 1650 km/s and c = 
1. Note also that only the wind has been allowed to be 
variable, i.e., we considered perturbations only outside the 
sonic point, located roughly at 20 km . 

In this way then, different and random amplitudes 
within the maximum range Sp/po G [—0.7, 0.7] are created 
for each individual shell. To allow for a temporal variabil- 
ity of the resulting profiles at "random" observation times 
(remember that our observations are typically separated 
by much more that one wind flow time) we performed, 
for each model, 10 simulations with different initialization 
of RAN and different locations of the contributing shells. 
The resulting 10 different profiles have been analyzed sub- 
sequently by means of the TVS method (an example is 
given in Fig, j^. 

Model series SS suffers from (at least) two major prob- 
lems. At first, the assumption of spherical shells might 
amplify the Ipv at all frequencies, since, especially in the 
wind lobes, there is only a weak chance that fluctuations 
will cancel out due to statistical effects. Second, our simu- 
lations define equal amplitudes of disturbance inside shells 
of equidistant range in velocity space. Accounting for the 
rather steep increase in velocity inside and the flat veloc- 
ity field outside, this means that the contributing volume 
per shell is strongly increasing with radius, which might 
give too much weight to disturbances in the outer wind. 
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Fig. 7. Left panel: "Observed" values of the numerator entering on (TVSintograi, asterisks) and of the fractional 
emission equivalent width (denominator of cn, crosses), compared with simulated quantities, as a function of mean 
wind density log < p >. The dotted curves correspond to the simulations with spherical shells and constant 6v spacing 
(SS1/SS2), and the other curves to the somewhat more realistic simulations accounting for broken shells and constant 
Sm, BS11/BS12 (fully drawn) and BS21/BS22 (dashed). Note that inside each series of simulations the maximum 
amplitude of density-fluctuations, 6p/po, is identical, i.e., independent on wind density. Right panel: As left, but for 
the fractional amplitude on (in %). The object numbers correspond to the entries given in Table |21 
In both panels, the special symbols correspond to the data resulting from our TVS-analysis of the 3-D models presented 
bv fHarriea (,2000,1 . cf. Sect. 16.21 (diamonds: spiral structure; triangles: clumpy structure). 
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Fig. 6. As Fig. ^ but for synthetic profiles correspond- 
ing to model BSll (broken shells, low amplitude) of 
HD 188 209, at M= 2.5 • lO-^Mg/yr (see text). 



To "cure" both problems, we have calculated four ad- 
ditional model series, which should be more realistic than 
the above ones. At first, the spherical symmetry is bro- 



ken by the following modification.^''' For the core-rays, we 
assume coherent shells (blobs), either of a relatively small 
lateral extent, Ap « i?*/10 (series BSl) or of a larger ex- 
tent, Ap — Ri, (series BS2). For each of the non-core rays 
(distributed roughly logarithmically), on the other hand, 
we assume different locations of the density variations per 
ray, to simulate the presence of broken shells. The latter 
modification results in a lower TVS particularly in the red 
part of the profiles, due to cancellation effects. Note that 
we have convinced ourselves that different distributions of 
non-core rays gave very similar results. 

In order to avoid the volume effect, instead of assuming 
5p/po as random, however constant per shell of thickness 
Sv — const, we now require that the random perturbations 
should occur in shells of equal mass, 

(^TOshcii = Airr'^pdr, (15) 

with roughly 50 (broken) shells per model. Inside each 
broken 6m shell, the density fluctuations are evaluated 
as above. For each of our simulations BSl and BS2, we 
have used two different values for the maximum ampli- 
tude, max(5p/po) = ±0.35 and ±0.7 (BS11/BS21 and 
BS12/BS22, respectively), which gives a fair consistency 
with the range of observed variability. 

Before we discuss the results in detail, let us already 
point out here the major outcome. Although the assump- 
tions inherent to the various model series (SS vs. BS) 

Remember that the radiative transfer is performed in the 
usual p — z geometry, with impact parameter p and height over 
equator z. The so-called core rays are defined by p < _R*, and 
the non-core rays passing both hemispheres of the wind lobes 
by p > Ri,. 



are rather dirterent, the results with respect to interest- 
ing quantities are fairly similar. The only difference con- 
cerns the distribution of the variability over the profile. 
For the spherical shells models, we find significant vari- 
ability on the red side, whereas for the broken shell model 
the variability extends to larger blue velocities, due to the 
increased influence of the shells in front of the disk (can- 
cellation effects in the lobes, see Fig.lHJ. 

In Fig. 13 we now compare the outcome of all our 
simulations with the observations ("conservative" values). 
On the left panel, we have plotted the numerator enter- 
ing On (lower set of curves), and the fractional emission 
equivalent width (FEEW, denominator of on, upper set 
of curves) as a function of mean wind density < p >. 
Obviously, series SS and BS give similar results, as al- 
ready noted. In particular, the results of models SSI and 
models BSll (lower dotted and fully drawn curves) are al- 
most identical, which shows that a large number (^ 160) 
of spherical shells (model SSI) can simulate the outcome 
of a model with a lower number (~ 50) of broken shells 
and a lower density contrast. Moreover, it seems that the 
"volume effect" discussed above is insignificant, simply be- 
cause Ha forms in the lower wind region. This similarity 
in the results points to a rather large probability that our 
results are robust and independent of the specific assump- 
tions. 

Interestingly, both the observations (except for the two 
objects with highest wind-density and more localized TVS, 
HD 190429A and HD 16 691) and all simulations roughly 
follow a power-law for both quantities, 

l0g(rT/S'intcgral) « fl + & log < p > 

log(FEEW) « c + d\og< p>, 

which immediately shows that our hypothesis of both 
quantities being proportional to each other fails. Note at 
first that the logarithmic dependence of the FEEW on 
log < p > can be readily understood if one remembers 
that the total emission equivalent width of Ha scales as a 
power-law of mean wind-density (cf. Puis et al. 1996), and 
that the integration range [vb , Vr] , entering the fractional 
equivalent width, is only weakly increasing with < p >, 
if scaled to v^o and evaluated on a logarithmic scale (see 
Fig.EJ. Since the TVS and its integral, on the other side, 
is also related to the mean wind-density (at least if the 
disturbances do not totally decouple from this quantity), 
the power-law dependence of this quantity on < p > can 
also be understood. The different and lower slope can be 
attributed to optical depth effects and the cancellation 
of fluctuations in the emission lobes (vs. the contribu- 
tions from core-rays), at least in our simulations (where 
we "know" the origin of the variability). In total then, on 
becomes a decreasing function of log < p >, 

log On « (a — c) -I- (6 - d) log < p >, b < d, (16) 

and since all our model series predict the same depen- 
dency, it is rather likely that this effect should be present 
also in more realistic simulations. In conclusion, we pre- 
dict that On becomes a decreasing function of mean wind 



density, even ij the disturbances (more precisely, their rel- 
ative amplitudes) are independent of < p >. The vertical 
off-set of this relation, on the other hand, depends strongly 
on the density contrast, i.e., on max((5p/po). All this sim- 
ulations, of course, refer to the case of fluctuations which 
are "globally" present, and will not explain effects from 
localized macro-structures such as CIRs. 

The actual and predicted behaviour of on is shown in 
Fig. d right panel. For lower wind densities, the slopes of 
the relations for numerator and denominator are rather 
similar (optically thin winds, wind emission dominated 
by core-rays), so that the predicted amplitude on re- 
mains roughly constant or is even increasing, whereas 
from log < p >« —14.0 the predicted decrease is obvi- 
ous. By comparison with observations, we find that almost 
all stars lie witin the range suggested by BS11/BS12 and 
BS21/BS22, i.e., correspond to differences in (relative) 
amplitude within a factor of two. There are only two out- 
liers, HD 210809 (#10) and particularly HD 18 409 (#15, 
at On ~ 50). Remember that the former star has been 
designated as a "peculiar" object in Sect. 15.11 whereas 
the strong deviation of HD 18 409 is more l ikely due to 
uncertainties in wind parameters ijRepolust et al. 2 00^ : 
the large value of cn is not due to a large TVS-integral, 
but due to a rather small value of its FEEW (located at 
log < p >w -13.64 and FEEW = 39 km/s in the left 
panel. 

Even for the two objects with the largest wind den- 
sities, which have been found not to follow the individ- 
ual relations for the TVS-integral and the FEEW (due to 
their rather localized variability), the results for on are 
consistent with the predictions. In our interpretation, this 
would mean that both stars have the same degree of activ- 
ity as the other stars, only in different and more localized 
regions. 

In summary, there are no indications of a dependence 
of wind activity on wind density, at least on basis of our 
present simulations; in our interpretation, the decrease of 
On is an artefact of the normalization, which (unfortu- 
nately) does not follow the same slope as the TVS-integral. 

Furthermore, from our simulations, we may also con- 
clude that the underlying disturbances (at given wind den- 
sity) may introduce a scatter up to a factor of two in the 
maximum amplitude. Of course, more realistic simulations 
are needed before a final statement concerning this point 
can be given. 

Further insight into the origin of the variability might 
be found from a comparison of observed vs. predicted ve- 
locity limits, Vb, Vr and particularly of their asymmetry, cf. 
Fig-IHl In the left panel of this figure, we compare these ve- 
locity limits (normalized to v^c , to obtain a unique scale) 
with our predictions, here with results from BS21 (bro- 
ken shells, low amplitude). Obviously, for stars with low 
and intermediate wind strength ([log < p >] < —13.2) 
and except for HD 210 809 (at log < p >= -13.6, with 
strongly asymmetric velocity limits) , our models do fairly 
reproduce the observed amount of increase of AV as a 
function of < p >. Since at largest wind densities we 
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Fig. 8. Left panel: Observed blue and red velocity limits, Vh,Vr (conservative and non-conservative values) in units of 
Woo , as a function of mean wind-density, compared with results from simulation BS21 (dashed). 

Right panel: Observed asymmetry, {vr — \vb\)l{vr + \vb\) (conservative values, black dots), compared with simulations 
(asterisks: BS21, squares: BSll, crosses: SS3). The object numbers correspond to the entries given in Tabled The 
diamond and triangle refer to the data resulting from our TVS-analysis of the 3-D spiral and clumpy model presented 
bv .Harries (,200a) . cf. Sect. lO 



have only two objects in our sample, it is not clear at 
present whether their discrepant behaviour is peculiar or 
not. Thus, except for the outliers, it might be concluded 
that the observed variability results from effects which are 
present everywhere in the wind, in accordance with our 
models. This conclusion seems to be also supported by 
the fact that HD 210 809 deviates from our predictions: 
for this star our observations have suggested the presence 
of large-scale wind disturbances which are localized rather 
than uniformly distributed over the wind volume. 

The right panel of Fig. |H1 however, shows also the 
shortcomings of our models. Plotted is the asymmetry of 
Wfc , Vr by means of the expression (vr — | | ) / (f r + 1 Wfc | ) (neg- 
ative values correspond to blue-to-red, positive values to 
red-to-blue asymmetries, respectively). Let us firstly dis- 
cuss the "theoretical" predictions. We have plotted the re- 
sults for model series BS21 (asterisks), BSll (squares) and 
for a model which has been constructed additionally for 
this comparison, namely a model with spherical shells and 
6m = const (SS3, crosses). All three models have the same 
maximum amplitude, max((5p//9o) — ±0.35. Not surpris- 
ingly, the last model shows the least asymmetry, whereas 
model BS21 shows the largest one, due to the large lateral 
extent of the assumed broken shells in front of the stellar 
disk. Note that most models show a blue-to-red asymme- 
try, even those with spherical shells, whereas in no case a 
redwards asymmetry is found. This predicition, of course, 
results from cancellation effects in the lobes, compared 
to the situation for core-rays. For large wind-densities, 
all models converge to small or even negligible asymme- 
try, because of the increasing influence of the lobes. For 
model SS3 (spherical shells), symmetry is reached earliest, 
roughly at log < p >— —13.5. 

With respect to the observations, the situation is as 
follows. At low and intermediate wind densities, only four 



objects display a significant asymmetry, mostly to the red, 
where the highest degree is found for HD 210809 (#10, 
see above). The majority of stars, however, show either 
a small degree of asymmetry (both to the blue and to 
the red) or behave symmetrically. The predictions of se- 
ries SS3 are closest to this behaviour. Let us point out 
already here that a symmetric TVS at low wind densities 
can be alternatively explained by models with co-rotating 
structures (e.g., CIRs), as will be discussed later on. For 
the two stars with larger wind densities, the predictions 
definitely deviate from what has been observed. 

In Sect. 15.3.21 it turned out that the mean amplitude of 
deviations. An, seems to be uncorrelated with mean wind 
density, at least for our sample. We speculated that if the 
sources of wind variability were (uniformly) distributed 
over the contributing wind volume, one should see a pos- 
itive correlation. By means of our simple models, we can 
check this conjecture now and may find additional con- 
straints on the origin of wind variability. 

Fig. ini shows the behaviour of the simulated mean 
amplitude as a function of mean wind density, for se- 
ries BS11/BS12, series BS21/BS22 and SS3, respectively. 
Actually, is strongly correlated with < p >, where the 
vertical offset is a function of (relative) amplitude Sp/po. 
On a first glance, we might conclude that the observed 
values are again (i.e., as we have found for on) consistent 
with the models, if we allow for a variation in amplitude 
of roughly a factor of two. 

A closer inspection, on the other hand, reveals that 
at least for two low density objects (HD 188 209, #11 
and HD 218 915, #14) there is the following problem. The 
mean amplitudes of these objects are consistent with our 
simulations with larger Sp/ po (BS22). Concerning the frac- 
tional amplitudes, on, however, they are consistent with 




Fig. 10. As Fig. n however for the time-series of synthetic Ha profiles from ^ Pup as calculated by iHarriesI (1200(1) . 
using two different 3-D models for the wind morphology. Left panel: wind with co-rotating, one-armed spiral density 
enhancement; right panel: clumped wind. Note the different distributions of the corresponding TVS. 
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Fig. 9. "Observed" values of the mean amplitude of devi- 
ations, An (conservative values) compared with simulated 
quantities, as a function of mean wind density log < p >. 
The bold curves correspond to series BS11/BS12, the 
dashed curves to series BS21/BS22 and the dotted curve 
to series SS3 (spherical shells, low amplitude), respec- 
tively. Each series shows a positive correlation with mean 
wind density. The object numbers correspond to the en- 
tries given in Tabled 



our simulations for lower Sp/pa (BS21, cf. Fig. Cj), which 



in turn produce a much too large bluewards asymmetry 
in the TVS (Fig. Remarkably, one of these objects 
(HD 218915, #14) even suffers from an observed red-to- 
blue asymmetry. 

Thus, from the comparison of mean amplitudes and asym- 
metry in the velocity limits we find a number of indications 
that at least two (from 3) low and the high density ob- 
jects cannot be explained by our simple models consisting 
of density fluctuations distributed everywhere in the wind. 
If we return to our original TVS analysis (Figs. |31 and 
the primary sources for these inconsistency can be clearly 
identified: (i) the strong variability of a central reversal 
at zero rest velocity, which gives rise to a rather large 
TVS within a narrow, symmetric velocity range for stars 
of weaker winds and (ii) the rather small extent of the ob- 
served variability, preferentially bluewards of the emission 
peak for stars of stronger winds. 

In our models, a large variation at rest wavelength can- 
not be reproduced, at least if we allow for fluctuations of 
both positive and negative amplitude. Such a behaviour 
might be explained by radially extended, coherent struc- 
tures in front of the disk, which would mimic a certain 
global increase of mass-loss by bringing the innermost part 
of Ha into emission. That our models can never repro- 
duce a (strongly) localized variability has been discussed 
already above. 



b.z. jU model simulations 

Recently. iHarriesI l|200(tl has published results for 3-D line- 
profile simulations of Ha for Pup, performed by means of 
his Monte Carlo stellar wind radiative transfer code. Two 
distinct models for the wind morphology have been consid- 
ered: one with a co-rotating (one-armed) spiral structure 
and another one consisting of a clumpy wind. In the first 
case the author examined the effect of one spiral stream- 
line of enhanced density on Ha, while in the latter one 
he considered random blobs propagating throughout the 
wind. 

In order to obtain an impression to what extent the 
outcome of our model analysis might be influenced by the 
fact that we consider 1-D instead of more realistic 3-D 
models, we analyzed the two different sets of synthetic 
profiles derived by Harries (kindly made available to us 
by the author), in the same way as for the time-series of 
our sample stars and for our own simulations. The corre- 
sponding results are shown in Fig. IIUI 

Apart from the impressive sequence of synthetic pro- 
files which allows to easily follow any evolution of wind 
structure in time a number of interesting features are no- 
ticeable. 

In the first place, the distributions of amplitudes for 
the two models are rather different: in the one-armed spi- 
ral model, the derived TVS is double-peaked and sym- 
metric with respect to the stellar rest frame, while in the 
clumped model it is single peaked with maximum am- 
plitudes concentrated on the blue. In addition, for the 
clumped model the velocity range of significant vari- 
ability shows a clear blue-to-red asymmetry (wb=-1670, 
Ur=1440 km s~^ ) while for the spiral model it is almost 
symmetrical with respect to the rest wavelength {vb=- 
2010, w^=2050 km s^^ ). 

Interestingly, also the parameters derived from the 
TVS analysis of the spectra for the 3-D clumped model 
are quite similar to the ones we have derived in terms of 
our 1-D broken shell (i.e., clumpy) models (see Figures 
andlSI), whereas the small differences might explained by 
the fact that the spectra from the 3-D simulations are free 
of noise, in contrast to the spectra from our 1-D models. 

These findings suggest that at least in the case of a 
clumpy wind the TVS signatures do not or little depend 
on the specific model, and that the results of our simple 
1-D simulations are comparable to those from more real- 
istic simulations. A final conclusion concerning this point 
can be drawn only, of course, when a larger set of 3-D sim- 
ulations for a variety of wind parameters, will performed. 

From the point of view of the observations, on the other 
hand, the above finding means also that the probability 
to obtain a symmetric TVS of Ha will be very low if the 
wind is clumpy, even if using snapshot observations. 



f. bpectral variability and the vvind Momentum 
Luminosity Relationship 

The most extreme case of snapshot observations occurs 
when performing a spectral analysis, since usually only 
one spectrum (or two consecutive ones) is/are available. 
Because of the spectral variability inherent to the objects 
under discussion, these individual data might be not very 
representative. In the following, we shall investigate this 
point with respect to derived mass-loss rates, with special 
emphasis on the question in how far the Wind Momentum 
Luminosity Relationship (WLR) of Galactic 0-type stars 
is infiuenced. 

To this end we assume (i) that the observed variabil- 
ity of Ha is only due to processes in the wind (i.e., the 
contribution of absorption Ipv is neglected) and (ii) that 
this variability is interpreted in terms of a variable mass- 
loss rate. Note that the latter assumption is inherent to a 
typical spectrum analysis, since due to the scarcity of the 
available data-set(s) a particular mass-loss rate is derived, 
which of course might be not representative. 

Following this philosophy, we derived upper and lower 
limits of the mass-loss rate for all our sample stars, from 
the two most extreme spectra present in the time-series, 
and calculated the corresponding limits for the modified 
wind-momentum rate. Due to the adopted simplifications 
(standard model, no clumping) the obtained values of 
M might be overestimated. 

The mass-loss rate estimates listed in Table El show 
that the observed variations in M range from ±4% of the 
mean value, for stars with stronger winds, to ±16% for 
stars with weaker winds. Since the accuracy of our de- 
terminations also depends on the strength of the wind 

- ±20%, for stars with Ha in emission, an d ±30%, for 
stars with pure absorption profiles (Markov a et alJl2004l) 

- we conclude that the observed variability of M does not 
exceed the errors of individual determinations and thus 
is not significant. This result seems somewhat astonish- 
ing, especially in those cases when drastic changes in the 
Ha profile shape have been observed. Note, however, that 
for not too low wind densities small changes in M give rise 
to large changes both in the profile shape and in the EW 
dPuls et al.lll99l . 

The derived amplitudes of the M variability have been 
used to cast constraints on the variability of the cor- 
responding wind momentum rates. The results indicate 
that wind variability affects the momentum rates by less 
than 0.16 dex, which is smaller than the error of individ- 
ual estimates (0.30 dex, Markova et al. 2004) and thus 
is insignificant again. In addition, and as it can be seen 
from Figure im the uncertainty caused by wind variability 
does not provide any clue to resolve the problem concern- 
ing the WLR of G alactic O-type stars being a function 
of luminositv class ijPuis et allbooilRepolust et al.ll2004t 
iMarkova et aPbOO^ ^ 

Finally, we would like to note that the assumption of a 
homogeneous and spherically symmetric wind (underlying 
our analysis) is in some obvious contrast with the presence 
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Fig. 11. WLR for our sample of Galactic 0-type stars. 
Error bars display the influence of Ha Ipv (interpreted in 
terms of M ) on the modified wind-momenta. 

of large-scale perturbations in the winds of some of our tar- 
gets, which have been suggested by different investigators 
on various occasions (see Sect.QJ. Consequently, problems 
may occur when trying to fit the Ha profiles with model 
calculations. Exemplary is the case of HD 210839, where 
we failed to obtain a good fit for a number of profiles of 
the time series. 

8. Summary, discussion and conclusions 

Although line-profile variability of Ha in 0-type super- 
giants has been known for a relativ ely long time (e.g., 
lRosendha]|[l973t IConti fc Fros3ll977|) . our survey is the 
first where this variability is detected and quantified in 
an objective and statistically rigorous manner for a large 
sample of stars. By means of a comparative analysis we 
furthermore were able to put a number of constraints on 
the properties of this variability as a function of stellar 
and wind parameters. 

The main results of our survey can be summarized as 
follows: 

To study the wind variability in O supergiants, as 
traced by Ha, we employed the TVS analysis, originally 
developed by iFuUcrton. Gies & Bolton (1996), however 
modified in such a way as to account for the effects of 
wind emission of various strength on the observed pro- 
files. As already predicted by these authors, the so-called 
fractional amplitude of deviations, which serves as a mea- 
sure for absorption Ipv, cannot serve as an indicator for 
the strength of wind variability. Instead, we introduced 
and used a new TVS parameter, cn, which measures the 
contribution of unit wind emission to the variability de- 
tected by the corresponding TVS. 

By means of this quantity we found that the observed 
variability of Ha per unit wind emission is a moderately 
decreasing function of mean wind density. This result 
might be interpreted as an indication that stronger winds 
are less active than weaker ones. However, at least on the 



basis of our present, simplistic Ime-pronlc simulations, this 
hypothesis cannot be supported. 

All sample stars show evidence of significant line- 
profile variations in Ha with mean amplitudes An be- 
tween 1 and 4% of the continuum strength. Since ab- 
sorption Ipv seems to be common among 0-type stars 
l)Fiillerton. Gies fc Boltonl[l99(ii) . these amplitudes can in 
principle be due to a combination of variations gener- 
ated both in the photosphere and in the wind. However, 
our analysis indicate that even at lower wind densities 
(log < p >« —14) contributions from wind effects have 
to be accounted for, and that for a number of stars wind 
variability is actually the dominating source. 

In their original study, FuUerton et al. concluded that 
the mean amplitude of deviations is an inappropriate 
quantity to compare Ipv in different lines of various stars 
because it does not account for the strength of the under- 
lying feature. Although true for the case of absorption Ipv, 
this conclusion has to be somewhat modified with respect 
to our investigations. By means of simple line- profile sim- 
ulations, the mean amplitude of deviations in Ha (defined 
by Eq. I10() has been predicted to be an increasing func- 
tion of mean wind density, at least in those cases where 
the disturbances in density were present everywhere in 
the wind. Interestingly, this prediction is not supported 
by our observations which gave no evidence of any cor- 
relation between mean amplitude and density. Moreover, 
it turned out that a comparison of observed mean and 
fractional amplitude with corresponding simulated quan- 
tities give different information, because of the different 
infiuence of either the contributing velocity range of sig- 
nificant variability or the contributing net wind emission, 
respectively. In so far, both quantities (i.e., fractional and 
mean amplitude) deserve their own right and have to be 
reproduced simultaneously by future models which claim 
to explain the observed variability. 

iFullerton. Gies fc BoltonI l)l996(l have found that the 
distribution of variability within an absorption profile can 
provide information about the cause of the variations, e.g., 
radial or non-radial pulsations. By means of a series of 
models with different wind morphology we showed that 
the TVS analysis of Ha can also give some insight into 
the structure of the wind, at least of its lower part where 
this line forms. Both models with spherically symmetric 
and with broken shells produce an Ha TVS with a blue- 
to-red asymmetry, if the structures are distributed uni- 
formly over the contributing wind volume and as long as 
the wind-density is not too large. The only difference in 
the outcome of the two kinds of models is the actual de- 
gree of the asymmetry: at the same mean density, shell 
models produce less asymmetry than those with broken 
shells (i.e., blobs). Moreover, this result does not seem to 
depend on whether snapshots or systematic simulations 
have been used. Note, however, that an asymmetric large- 
scale, long-lived co-rotating structure will always produce 
a symmetric TVS z/the whole period of rotation is covered 
by observations. 



A comparison between the observations and line- 
profile simulations (assuming variations caused by coher- 
ent shells/blobs distributed everywhere in the wind) re- 
vealed that for stars of intermediate wind density the pa- 
rameters of the Ha TVS are consistent with those from 
the models if one allows for a scatter within a factor of two 
for the maximum (relative) amplitude of the disturbances 
at given wind density. 

On the other hand, the established disagreement at 
lower and higher wind densities might simply point to the 
presence of wind effects which are not accounted for in 
our simulations. For example, to reproduce the variations 
observed in Ha at lower wind densities one might suggest 
the presence of a radially extended, coherent structure in 
front of the disk which can mimic variations in the global 
mass-loss rate. The major problem concerning the (two!) 
stars with stronger winds relates to the fact that the ob- 
servations yield a TVS which is localized preferentially 
bluewards of the rest frame (inside a rather narrow ve- 
locity range), whereas our simulations always produce an 
almost symmetric TVS instead. This discrepancy might 
be again interpreted in terms of a rather confined region 
of variability close to the star. 

Our analysis shows that significant variations in Ha oc- 
cur below 0.3 V rp ■ This estimate is a bit higher than the 
value derived bv lKaper et al.l^|l997^ . namely 0.2 . Since 
the velocity extent of the observed variations depends 
strongly on < p > and since our sample includes a couple 
of stars with strong winds (Ha completely in emission) 
which are missing in the sample of Kaper et al., such a 
discrepancy is quite natural. Moreover, the velocity limits 
as derived by us, converted to units of physical space, cor- 
respond to 1.4 to 1.5 Ri, , in good agreement with results 
from theoretical calculations with respect to the outer lim- 
its of Ha line-formation in O supergiant winds. 

Interpreted in terms of variable mass-loss rates, the ob- 
served (partly extreme) variations in the Ha wind emis- 
sion indicate variations in Mof ±4% of the mean value 
for stars with stronger winds and of ± 16% for stars with 
weaker winds. The ratio of maximum to minimum mass- 
loss rate determined over the time interval present for our 
sample ranges from 1.08 to 1.47, with a tendency that 
weaker winds show larger values. The mean ratio aver- 
aged over the whole samp le is 1.22, in good agree ment 
with the value provided bv lPrinia fc HowarthI l)l986l) . 

The consequences of wind variability with respect to 
the wind momentum rate of the sample stars is smaller 
than 0.16 dex and hence not significant, taken the indi- 
vidual errors inherent to any M deter mination. This re- 
sult agrees well with an investigation bv lKudritzkil l|l999l) 
who reported 0.15 dex as the error in the wind momen- 
tum rate of one(!) A-supergiant, HD 92 207, introduced by 
wind variability. Thus, we conclude that wind variability 
in O supergiants does not seem to affect the main con- 
cept of the WLR although it might contribute to the local 
scatter by moving individual points (up and down) with 
respect to each other. 



b maily, it may be worth noting that tlic Ha prohles 
of the sample stars derived throughout our observations 
are quite similar, both in shape and strength, to those 
obtained by other investigators in vari ous epochs (e.g., 
B,osendhaI.197-3:.Conti Frost 1977: Scuderi et alJlliB 
lEbbetal Il982t lUnderhilll Il995l: iKaper et alJ Il997l) . This 
finding suggests that the atmospheres of our targets are 
not subject to long-term variability. On the other hand, 
no information about previous Ha observations concern- 
ing the stars BD-f 56 739 and HD 338 926 was found in the 
literature. 
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